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ABSTRACT 


Data  collected  for  calculating  electron  content  of  the  ionosphere 
by  the  Farediy  rotation  technique  can  be  used  to  calculate  the  magnetic 
dip  angle  at  ionospheric  heights  near  the  geomagnetic  dip  equator.  In 
this  report  the  magnetic  dip  angle  at  ionospheric  height  was  determined 
at  the  position  where  the  angle  between  the  ray  path  from  the  satellite 
to  the  ground  station  and  the  geomagnetic  field  vector  is  90  (trans¬ 
verse  position). 

The  magnetic  dip  angle  calculation  was  based  on  the  assumption  of 
a  constant  height  of  350  km  for  the  centroid  of  the  ionospheric  electron 
density  profile.  The  spherical  harmonic  analysis  is  used  to  verify  the 
experimental  results.  The  results  compare  well  with  the  surface  mag¬ 
netic  dip  angle,  measured  in  1956-1966  by  the  Geodesy  Department  of  the 
Ministry  of  Defense,  Thailand,  and  indicate  thet  the  magnetic  dip  angle 
at  ionospheric  height  is  very  close  to  the  surface  value. 

These  data  permit  estimation  of  the  position  of  the  geomagnetic 
dip  equator  at  the  ionospheric  height.  The  estimated  position  of  this 
equator  is  directly  above  the  geographic  latitude  of  9.30°N.  This 
latitude  intersects  the  southern  peninsula  of  Thailand  about  *80  km 
south  of  Bangkok  (near  Surat  Thani).  This  estimated  value  compares 
well  with  the  surface  geomagnetic  dip  equator  and  indicates  that  the 
dipole  field  is  a  good  model  for  Southeast  Asia  in  the  vicinity  of 
longitude  100°fc. 
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I  INTRODUCTION 


The  Electronic s  Laboratory  of  the  Military  Research  and  Development 
Center  (MRDC)  Bangkok,  Thailand  was  set  up  as  a  joint  Tha i 1 and-Un i ted  St ates 
organization  for  the  study  of  communication  problems  in  equatorial  areas. 

One  of  the  major  objectives  of  communication  research  in  equatorial  areas 
is  the  study  of  the  ionosphere.  To  increase  understanding  of  the  i onospher  ic 
characteristics,  records  of  the  Faraday  rotation  on  40-  and  41-MHz  signals 
from  the  S-66  (EXPLORER  22)  radio  beacon  satellite  have  been  accumulated  at 
the  Bangkok  MRDC  Electronics  Laboratory ,  1,4  and  at  four  remote  sites  located 
north  and  south  of  Bangkok.  The  projection  of  a  typical  satellite  path  on 
the  earth  over  these  sites  is  shown  in  Fig.  1.  Thus  data  for  sites  in  the 
rouge  of  geographic  latitude  between  about  7°N  and  1 ')  N  were  collected. 

Data  collected  for  calculating  electron  content  of  the  ionosphere  by 
the  Faraday  rotation  technique  can  be  used  to  calculate  the  magnetic  dip 
angle  above  the  earth'.-:  surface  near  the  magnetic  equatorial  region.  The. 
application  of  this  technique  will  help  in  locating  the  geomagnetic  dip 
equator,  which  passes  through  Thailand  to  the  south  of  Bangkok.  The  tech¬ 
nique  provides  information  foi  comparison  of  the  above-surface  magnetic  dip 
angle  with  the  surface  magnetic  dip  angle  available  from  measurements  taken 
at  the  earth's  surface  mi  Thailand,5  with  other  published  values,6,7  and  with 
the  theoretical  value  brsed  on  the  spherical  harmonic  analysis.  Thus  this 
work  will  aid  in  studying  radio  propagation  effects  related  to  geomagnetic 
control  o'  the  ionosphere  that  are  important  to  communication  in  Southeast 
Asia.  Fo  example,  consideration  of  the  magnetoion ic  theory  ied  llagn  to  ex¬ 
plore  the  optimum  orientation  for  linearly  polarized  antennas  used  on  sliotc 
ionospheric  paths  near  the  geomagnetic  dip  equator.8  The  present  work  will 
permit  a  better  estimate  of  where  to  use  linear  antenna  (e.g.,  dipole)  ori¬ 
entation  and  where  a  diversity  scheme  based  on  antenna  orientation  will  work 
(owing  to  the  effects  of  the  earth’s  magnetic  field  in  the  ionosphere). 

This  work  also  provides  information  for  studying  the  geomagnetic  controlled 
HF  spectrum,  and  that  can  be  used  to  determine  where  post-sunset  spread-F 
conditions  are  most  likely  to  occur  and  hence  cause  HF  communication  outage. 


*  Reference*  *re  litted  *i  tlie  end  of  the  report. 
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I  Ins  report  presents  the  method  of  calculations  and  results  ler 
selected  Faraday  rotation  data  from  satellite  posses  observed  by  the  re  - 
ceding  stations.  The  assumptions  involved  in  calculations  are  given  and 
the  possible  source,  of  error  are  discussed. 
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II  THEORETICAL  BACKGROUND 


When  an  electromagnetic  wave  propagates  through  the  ionosphere  from 
a  satellite,  its  resultant  polarization  progressively  rotates,  under  the 
influence  of  the  magnetic  field,9  and  this  phenomenon  is  called  Faraday 
rotation.  Figure  2  illustrates  the  concept  of  Faraday  rotation.  The 
satellite  is  shown  in  t h ree  pos i t i ons  (d  >  90  ,  0  -  90  ,  and  d  <  90°, 
where  r‘  is  the  angle  between  the  signal  ray  path  and  the  earth's  magnetic 
field),  and  the  ionosphere  between  the  satellite  and  the  ground  station 
is  indicated.  The  centroid  height  of  the  ionospheric  electron  density 
profile  is  shown  as  hc  .  At  one  time  during  the  passage  of  the  satellite, 
the  signal  ray  path  is  pe; pendi c u  I  a r  to  the  magnetic  field  (0  -  90°). 

•Sale  I  1  i  i  _*  transit  of  this  unique  point,  called  the  transverse  position 
TQ ,  usually  leaves  r.  visible  signature  on  the  Faraday  rota,  ion  fading 
records  for  stations  located  near  the  geomagnetic  equ'lor,  An  example 
of  sucl.  a  sigiature  is  shown  on  the  failing  record  of  lig.  3. 
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LOCAL  TIME 


FIG.  3  SAMPLE  OF  SATELLITE  FADING  RECORD 

By  using  t' ••  “stellite  ephemeris  data  for  satellite  position  informa¬ 
tion  ar d  the  observed  time  of  Tq,  it  is  possible  to  calculate  the  magnetic 
dip  angle  from  magnetic- field  and  ray-path  geometry  considentions.  The 
value  of  dip  angle  thus  obtained  is  considered  to  apply  at  a  certain  height 
in  the  ionosphere.  This  height  is  chosen  by  the  following  considerations: 
The  electron-density  distribution  profile  of  the  ionosphere  based  on  the 
Chapman-type  distribution  has  found  wide  use  in  modeling  the  topside 
ionosphere.  Although  several  very  questionable  assumptions  were  made  in 
the  derivation  of  Chapman's  formula,  the  general  shape  of  the  electron- 
dens'ty  distribution  with  altitude  as  given  by  Chapman’s  equation  still 
gives  a  reasonable  fit  of  the  actual  distribution  profile  as  we  know  it.® 
From  the  shape  of  the  ionospheric  electron-density  distribution  it  is 
assumed  that  most  of  the  Faraday  rotations  occur  near  the  peak  of  this 
distribution  in  the  F  layer.  9,11  This  study  assumes,  as  other  studies  have, 
that  the  main  effect  of  the  earth’s  magnetic  field  is  localized  near  the 
centroid  of  the  e  lect  ron- den  si  ty  profile.11,12  The  centroid  height  for  the 
area  of  interest  is  obtained  by  fitting  the  Faraday  rotation  electron  con¬ 
tent  data  to  a  Chapman  model  of  the  ionosphere.  So  the  centroid  height 
in  the  ionosphere  is  determined  to  be  at  a  level  approximately  50  km 
above  the  height  of  maximum  dens  i  ty,  9,11  which  can  be  determined  from  a 
true-height  analysis  of  vert i ca i - incidence  ionograms  obtained  near  the 
location  and  time  of  v.he  satellite  passages. 


6 


Ill  MAGNETIC  DIP  ANGLE  CALCULATION 


The  purpose  of  this  section  is  to  describe  the  method  of  calculation 
of  magnetic  dip  angle  at  ionospheric  heights  near  the  geomagnetic  equa¬ 
torial  region. 

The  Faraday  rotation  technique  used  to  determine  the  magnetic  dip 
angle  at  the  centroid  height  of  the  ionospheric  electron  density  profile 
was  applied  at  the  position  whete  the  TQ  (transverse  position)  occurred. 

As  illustrated  in  Fig.  3,  the  local  time*  was  recorded  on  the  chart  to 
relate  the  observed  data  to  the  physical  position  of  the  satellite  in  its 
orbit.  The  unique  point,  TQ,  is  plainly  visible  in  this  chart.  The 
satellite  position  and  height  at  the  time  corresponding  to  TQ  can  be  ob¬ 
tained  from  ephemeris  data  provided  by  the  U.S.  National  Aeronautics  and 
Space  Administration  (NASA). 

For  the  calculations  of  this  section  the  wave  frequency  was  assumed 
to  be  so  high  that  refraction  can  be  neglected  (that  is,  the  wave  frequency 
is  much  higher  than  the  maximum  plasma  and  collision  frequencies).  A 
profile  view  of  lie  satellite  and  the  preund  station  is  shown  in  Fig.  4. 

It  can  be  s-eii  that  the  coordinates  o*'  the  sub-satellite  point  (the  pro¬ 
jection  of  satellite  on  the  earth),  the  sub- i onospher i c  point  (the  pro¬ 
jection  on  the  earth  of  the  intersection  of  the  ray  path  between  satellite 
and  ground  station  with  the  centroid  of  the  ionospheric  electron  density 
profile)  and  ground  station  are  A.  ;  0 1 ,  \  ;  t?r,  \  ;  geographic  latitude 

and  longitude,  respectively.  The  sa?  rllite  height  is  h  .  The  earth’s 
radius  is  denoted  by  fl  and  the  ionospheric  centroid  height  is  h.c. 

At  the  position  where  the  signal  ray  path  from  the  satellite  is  per¬ 
pendicular  to  the  geomagnetic  field  (T  position),  the  magnetic  dip  angle, 

/ at  the  ionospheric  centroid  height  can  be  calculated  from  the  expres¬ 
sions  given  by  the  following: 

/  8  =  tan  1  (cot  cos  Dfl)l  (1) 

*  Local  lime  waft  established  using  a  Geneial  Hadio  11 1 5B  Frequency  Standard  adjusted  using  phase  referenced 
to  the  VLF  transmission  from  GBR  in  England  on  16  kHz.  Time  ticks  generated  by  a  GR  1I23A  Sync  Tonometer 
from  ihis  standard  were  compared  daily  with  time  transmissions  from  WWVH  ard  J.IY  on  10  and  15  IMt.  The 
resulting  time  estimate  is  accurate  to  +5  mi  1 1 isec onda . 

^  The  derivation  is  given  in  Appendix  A. 
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FIG.  4  GEOMAGNETIC  FIELD  —  RAY  PATH  PROFILE 
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The  magnetic  dip  angle  obtained  from  Kq.  (1)  t°»'  lll,‘  ionospheric 

centroid  height  corresponds  to  the  sub- i  oiiospher  i  c  point,  f ; ,  and  A,.  Tb 
sub-ionospheric  latitude,  l,\,  can  be  obtained  from  Eq.  (6),  and  the  sub- 
ionospheric  longitude,  A,,  can  be  calculated  by  the  following  equations: 
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(9) 


These  relationships  between  the  satellite  and  the  ground  station  at 
the  T0  position  can  be  used  to  calculate  the  magnetic  dip  angle  at  the 
ionospheric  centroid  height  from  the  experimentally  obtained  Faraday  ro¬ 
tation  data. 
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IV  SPHERICAL  HARMONIC  MODEL  CALCULATION 


The  geomagnetic  field  can  be  represented  theoretically  by  mathemat¬ 
ical  functions.  ',l3  The  spiericni  harmonic  expansion  is  chosen  to  represent 
a  theoretical  model  of  the  earth’s  magnetic  field,  since  its  functions  are 
suitable  for  the  mathematic  expression  of  arbitrary  distributions  on  the 
surface  of  any  sphere. 

The  spherical  harmonic  function  to  be  considered  here  to  serve  for 
the  earth’s  magnetic  field  representation  is  based  on  several  basic  assump¬ 
tions,  since  the  fundamental  causes  for  the  earth’s  field  have  not  been 
satisfactorily  explained.  Gauss assumed  that  there  is  no  magnetic  matter 
near  the  ground  and  no  electric  current  passing  from  the  atmosphere  to  the 
ground.  Mence,  if  we  assume  that  there  is  no  external  source,  all  the 
currents  producing  the  magnetic  field  are  inside  the  earth's  surface.  Then 
in  any  solution  to  find  the  J i s wr i but i on  of  magnetic  field  in  the  region 
that  contains  no  current,  the  curl  of  the  magnetic  field  vector,  B,  at  any 
point  must  (by  Maxwell's  equation)  be  zero: 

V  x  B  ,  0  .  (10) 

This  equation  implies  the  existence  of  the  gradient  of  the  magnetic 
potential  function,  <p,  as 

B  =  -V<p  .  (11) 

By  Maxwell’s  equation,  the  divergence  of  B  is  identically  zero: 

V  •  B  =  0  .  (12) 

Substituting  Eq.  (11)  into  Eq.  (12)  the  magnetic  potential  function 
•satisfi  .s  Laplace's  equation: 

72sfc  =  0  (13) 

The  solutions  of  this  equation  in  spherical  coordinates  are  known  as 
spherical  harmonic  functions.  13,14  Figure  5  shows  the  spherical  coDrdinate 
system,  colatitude  (t1)  increasing  southward,  longitude  (/V)  inc  reas  i  ng  east¬ 
ward,  and  radial  distance  (a)  increasing  outward.  The  magnetic  potential 
function  at  any  point  can  be  expressed  .is: 
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Here  Rf  denotes  the  radius  ol'  the  earth  fequal  to  6370  km),  4*  and 
fi"  are  the  numerical  harmonic  coefficients;  and  P*(cos  H)  is  ,..ie  Schmidt 
orthogonal  function  of  fegree  n  and  order  m  with  m  £  n  >_  1.  The  Schmidt 
function  can  be  expressed  as  numerical  multiples  of  the  associated 
Legendre  polynomial,  Pn  t(cos  0) ,  as 
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By  Eq.  (11)  the  magnetic  induction  can  be  written  by  the  differential 
equation 
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The  earth’s  magnetic  induction,  B,  at  any  point  is  a  vector  and  its 
modulus  is  called  total  intensity  or  total  field.  It  is  noted  that  the 
earth’s  magnetic  field  has  been  tacitly  treated  mathematically  as  if  it 
were  constant  in  time.  This  total  field  is  defined  usually  by  the 
Cartesian  components  as  positive  to  the  northward  component  (X),  eastward 
component  ()’),  and  downward  component  (Z) ,  as  shown  by  the  coordinates  of 
point  P  in  Fig.  5. 
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Hence,  the  Cartesian  components  of  total  field  can  be  obtained: 


1 

1 

a  sin  6  ^ 


f  17) 


The  total  field  can  be  obtained  by  B  a  (H2  *  Z2)^  where 
H  *  ( X 2  +  Y2)^.  The  magnetic  dip  angle  7 B  =  arc  tan  Z/H. 

The  spherical  ha  monic  function  wts  calculated  using  the  CDC-3100 
electronic  computer.  The  program  was  performed  for  the  first  six  harmonics 
using  the  spherical  harmonic  coefficients  available  to  the  author.  15 
Figure  6  shows  a  typical  plot  of  the  computed  total  field  decreasing  with 
the  increasing  of  the  higher  degree  of  n  at  Bangkok.  The  ratio  (fl(,/a)"+l 
*  l  Ec .  (14)  is  the  altitude  decay  factor.  This  factor  indicates  the  de- 
crease  of  total  field  with  the  increase  of  altitude;  the  greater  the  alti¬ 
tude  up  to  the  height  of  ionosphere,  a  =  Rt  *  h c,  the  higher  the  degree  of 
n  that  can  be  neglected. 

Values  of  magnetic  di  p  angle  computed  at  various  altitudes  over  the  re¬ 
ceiving  stations  are  shown  in  Table  i .  These  results  show  that  the  magnetic  dip 
angle  even  at  an  altitude  of  400  km  is  very  close  to  the  surface  values.  Figure  7 
presents  a  map  of  the  surface  magnetic  dip  angle  of  Thailand  prepared  from 
this  computer  program.  The  spherical  harmonic  equations  of  magnetic  field 
components  for  the  first  six  harmonics  are  shown  in  detail  in  Appendix  B. 


Table  I 

COMPUTED  MAGNETIC  DIP  ANGLE 


LOCATION 

MAGNETIC  DIP  ANGLE 

Longitude 

Lat  it  ude 

Surface 

Value 

BE 

BSD 

s  400  km 

Chiang  Mai 

98.97°E 

18.82°N 

22.32°N 

22. 13°N 

22. 10°N 

22 . 06  °N 

Nakhon  Sawan 

100. 18°E 

15.66  °N 

14.98°N 

14.88°N 

14.86°N 

14.84°N 

Bangkok 

100.57°E 

13.73°N 

10.31°N 

10. 30°N 

10.29°N 

10.28°N 

Prachuap 

99.80°E 

1 1 . 80  °N 

5.91°N 

5. 90  °N 

5.90°N 

5.90°N 

Songkhla 

100.62°E 

7.20°N 

4.96°S 

4.87°S 

4.86°S 

4.85°S 

TOTAL  FIELD 


V  RESULTS 


A.  Trunsvf  se  Position  Magnetic  Dip  Angle 

Tlie  magnetic  dip  angle  over  each  receiving  site  was  calculated  at 
the  position  where  T0  occurs.  The  magnetic  dip  angle  was  assumed  to 
apply  at  the  centroid  height  of  the  ionospheric  profile  C h  approximately 
350  km).  This  height  appears  to  be  realistic  from  a  limited  amount  of 
t rue  -  he i ght  analysis  performed  on  typical  daytime  and  nighttime  ionograrns 
at  Bangkok  . 16 

The  results  of  the  magnetic  dip  angle  calculations  from  Faraday 
rotation  observations  of  live  receiving  sites  are  shown  in  Tables  II 
through  VI.  The  mean  and  variance  values  are  summarized  in  Table  VII. 

The  r.ean  values  of  magnetic  dip  angle  are  plotted  in  Fig.  8.  This  figure 
represents  the  magnetic  dip  angle  at  the  altitude  of  350  km  over  the 
range  of  geographic  latitude  between  about  6°  to  21CN  near  the  100°E 
longitude. 


Table  II 


MAGNETIC  DIP  ANGLE  AT  ALTITUDE  350  km  OVER  GEOGRAPHIC 
LOCATION  20.331*  .  99.21°E  (tHIANG  MAI  SITE  DATA) 


date 

SATELI.ITE 

REVOLUTION 

NUMBER 

TRANSVERSE 
POSITION  <T0! 
LOCAL  TIME 

MAGNETIC 

DIP  ANGLE 

29  Nov  1966 

107 24A* 

2048  r 15 

25. 42°N 

2  Dec  1966 

H)7o5A 

2025: 12 

25. 38 °N 

8  Dec  1966 

108  47 A 

1939:21 

25:  46 °N 

15  Jan  1967 

1136 5A 

1413:46 

25.  53°N 

18  Jan  1%7 

11406A 

1350:53 

25.  54°N 

21  Jan  1967 

11 447 A 

1327:58 

25.  41  °N 

26  Jan  1967 

11508D* 

0043:05 

25. 53“N 

27  Jan  1967 

11529A 

1242:10 

25.  52°N 

The  letters  A  ind  D  indicate  sscending  and  descending  (north- 
snd  south-going)  sstel!  otsses,  respectively. 


DATE 

SATELLITE 

REVOLUTION 

NUMBER 

TRANSVERSE 
POSITION  (T0) 
LOCAL  TIME 

MAGNETIC 

DIP  ANGLE 

28  C.-t  1966 

1027  2 A 

0059:56 

16.67°N 

29  Oct  1966 

10299D 

1237:12 

16 . 49°N 

31  Oct  1966 

1031 3A 

0036:57 

16.62°N 

3  Nov  1966 

1035 4A 

0013:57 

16. 43°N 

4  Nov  1966 

1037  4D 

1151:23 

16 . 5b°N 

5  Nov  1966 

10395A 

2350:57 

lo.64“N 

8  Nov  1966 

10436A 

2  327  :  56 

10 . 66°N 

11  Nov  1966 

10477A 

2304:57 

16.05°N 

Table  IV 

MAGNETIC  DIP  ANGLE  AT  ALTITUDE  350  km  OVER  GEOGRAPHIC 
LOCATION  14. 38°N,  I00.60°E  (BANGKOK  SITE  DATA) 


DATE 

SATELLITE 

REVOLUTION 

NUMBER 

TRANSVERSE 
POSITION  <rn) 
LOCAL  TIME 

MAGNETIC 
DIP  ANGLE 

25  Dec  1964 

105IA 

2126:23 

12.02°N 

10  Feb  1965 

I693A 

1456:42 

12.70°N 

18  Sept  1965 

47I2A 

0842:40 

12.07°N 

27  Oct  1965 

525ID 

1459:22 

12. 25°N 

30  Oct  1965 

5292D 

1436:35 

12. 24°N 

5  Nov  1965 

53740 

1351:11 

12. 05°N 

8  Nov  1965 

5415D 

1328:13 

12. 13°N 

10  Nov  1965 

5  4  36A 

0125:26 

12.67°N 

13  Nov  1965 

54"  7  A 

0102:43 

12.64°N 

16  Nov  1965 

5518A 

0039:59 

12.771S 

22  Nov  1965 

5600A 

1654:31 

12. 68°N 

1 3  J  an  1 966 

6324A 

1637  :  01 

12. 00°N 

16  Jan  1966 

6365A 

1614:23 

12. 11°N 

16  Mar  1966 

7 171 A 

0809:49 

12.24°N 

Table  V 


MAGNETIC  DIP  ANGLE  AT  ALTITUDE  350  km  OVEH  GEOGRAPHIC 
LOCATION  12.  14°N,  100.47°K  (PRACIH  AP  SITE  DATA  > 


DATE 

SATELLITE 

REVOLUTION 

NUMBER 

TRANSVERSE 
POSITION  <r0> 
LOCAL  TIME 

MAGNETIC 
DIP  ANGLE 

4  June  1%0 

827  0D 

0838:18 

6. 20°N 

7  June  1%(> 

831  ID 

0815:33 

6 . 1 7  "N 

9  June  1%6 

8  3  46  A 

2039:35 

6 . 20  °N 

12  June  1%() 

8  387 A 

2016:39 

6.50°N 

19  June  l%b 

84751) 

0644:  12 

5. 95°N 

20  June  1 9t)h 

84%A 

1840  56 

b .  37  “N 

22  June  1%0 

85161) 

0621:22 

6.00°N 

23  June  1%6 

8537A 

1817:35 

6. 21°N 

2b  June  1%6 

8  578  A 

1754:41 

6. 19“N 

29  June  l%b 

86  19A 

1731:44 

6 . 17  °N 

Tabe  VI 

MAGNETIC  DIP  ANGLE  AT  ALTITUDE  350  km  OVER  GEOGRAPHIC 
LOCATION  6 . 75°N,  100.66 °E  (SONGXHLA  SITE  DATA) 


DATE 

SATELLITE 

REVOLUTION 

NUMBER 

TRANSVERSE 
POSH  ION  (r„> 
LOCAL  TIME 

MAGNETIC 
DIP  ANGLE 

4  No  ■  1966 

1037  4D 

1154:56 

5.83°S 

6  Nov  19b6 

10395A 

2347:39 

5 . 97  °S 

8  Nov  1966 

10436  A 

2324:40 

6.00rS 

11  Nov  1966 

10477  A 

2301:38 

6.00°S 

14  Nov  1966 

1 05 18 A 

2238:41 

6. 02°S 

17  Nov  1966 

10559A 

2215:37 

6.00°S 

20  Nov  1966 

10600A 

2 1 52 :  37 

5. 99°S 

3  Dec  1966 

1077  ID 

0738  :  37 

5.82°S 

16  Dec  1966 

10956A 

1758  :  21 

5. 92°S 

19  Dec  1966 

10997 A 

1735:21 

5. 86°S 
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Table  VII 

MAGNETIC  DIP  ANGLE  AT  ALTITUDE  350  km 


SATELLITE 

RECEIVING 

STATIONS 

SUB- IONOSPHERIC  POINT 

MEAN 

MAGNETIC 

DIP 

ANGLE 

<v 

VARIANCE 

Longi tude 

<Ai ) 

L*t i t ude 

<*  > 

Chiang  Mai 

99. 21°E 

20. 33°N 

25. 47°N 

O.Ofc 

Nakhon  Sawan 

100 ,60°E 

16. 54DN 

16.59°N 

0.08 

Bangkok 

1 00 . 0  0°E 

1 4. 38UN 

12. 33°N 

0.08 

Prachuap 

100. 47°E 

12. 14°N 

6 . 20°N 

0.02 

Songkhla 

100 .60°E 

6.75”N 

5. 94°S 

0.07 

The  principal  feature  of  magnetic  dip  angle  is  its  linear  variation 
with  the  geographic  latitude  near  the  geomagnetic  dip  equator.  The 
estima.  _d  zero  magnetic  dip  occurred  directly  above  the  geographic- 
latitude  of  approximately  9.30  N  (about  480  km  south  of  Bangkok). 

B.  Comparison  of  Experimental  and  Theoretical  Determination  of 
Magnetic  Dip  Angle 

To  check  the  results  of  the  magnetic  dip  angle  at  the  ionospheric 
height  obtained  from  the  experimental  Faraday  rotation  records,  model 
calculations  were  made.  The  spherical  harmonic  function  was  chosen  to 
represent  -.he  theoretical  model  of  the  geomagnetic  field  distribution. 

In  addition,  the  surface  magnetic  dip  angle  measured  in  1956-1966,  by 
the  Geodesy  Department  Ministry  of  Defense,  Thailand5  and  the  published 
NOO-CGS  1965  Chart  1 7 00 6  were  used  to  verify  the  applicability  of  the 
results. 

The  accuracy  of  the  spherical  harmonic  function  in  modeling  the 
geomagnetic  field  distribution  depends  primarily  on  the  number  of  har¬ 
monics  used  The  significant  feature  of  the  geomagnetic  field  can  be 
accurately  obtained  by  involving  only  the  first  six  harmonics.  However, 
this  iib  t  hema  t  i  ca  1  function  cannot  give  any  representation  of  the  actual 
local  variation  of  the  geomagnetic  field. 
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The  comparison  of  the  magnetic  dip  angle  obtained  from  Faraday  rota 
tion  technique  with  those  derived  from  spherical  harmonic  function  and 
other  surface  vc'  es  is  shown  in  Table  VIII.  The  mean  Faraday  rotation 
technique  values  of  tiie  magnetic  dip  angle  were  taken  from  the  data  of 
Table  VII.  The  surface  measured  values  were  obtained  directly  or  inter¬ 
polated  for  the  corresponding  sub- i onospher i c  points  in  Table  VII.  The 
surface  published  values  were  interpolated  from  NOO-CGS  1965  Chart  1700. 
The  general  feature  of  Table  VIII  is  that  the  magnetic  dip  angle  at  the 
ionospheric  height  from  Faraday  rotation  technique  agrees  surprisingly 
well  with  those  surface  values. 


Table  VI 1 1 

COMPARISON  OF  MAGNETIC.  DIP  ANGLE 


GEOGRAPHIC 

LOCATION 

FARADAY 

ROTATION 

TECHNIQUE 

(A  ~  350  km) 
e 

SPHERICAL 

HARMONIC 

Uf  =  550  km) 

MOD 

SURFACE 

MEASURE¬ 

MENTS 

NOO-CGSj 

1965  CHART 
SURFACE 
VALUE 

Latitude 

Longitude 

20.33°N 

99.21°E 

25.47°N 

25.29°N 

25.80°N 

25.5°N 

16. 58°N 

100.60°E 

16.  59°N 

17. 19°N 

lb.80°N 

17 . 00°N 

14. 38 °N 

100.60°E 

12.33°N 

12.07°N 

12. 00°N 

12. 00°N 

12.14°N 

100. 47°E 

6.20°\ 

6.8°N 

6 . 50°N 

6. 5°N 

6.75°N 

100.66°E 

5.94°S 

5.  96°S 

5 . 90°S 

6.0°S 

Department  of  Geodesy,  Ministry  of  Defense,  Thailand. 


C.  Geomagnetic  Dip  Equator 

A  further  result  of  this  study  is  that  the  data  are  analyzed  to 
obtain  the  location  of  geomagnetic  din  equator  at  the  ionospheric  height. 
In  this  report  the  term  geomagnetic  dip  equator  means  the  locus  where  the 
magnetic  dip  angle  is  zero  Table  IX  shows  the  estimated  position  of 
geomagnetic  dip  equator  compared  with  estimates  by  other  workers.  All 
geographic  locations  of  the  geomagnetic  dip  equator  contained  in  Table  IX 
were  obtained  directly  or  interpolated  from  the  data  available  to  the 
author.  In  this  table  the  geomagnetic  dip  equator  estimates  agree  well, 
and  the  location  of  this  equator  at  the  ionospheric  height  is  very  close 
to  the  surface  value. 


Table  IX 


ESTIMATED  POSITION  OF  GEOMAGNETIC 

DIP  equator  near  ioo°e 


— 

SOURCE 

ALTITUDE 
(  km) 

GEOGRAPHIC 

LATITUDE 

Author  (1965-67) 

350 

9.30N 

TEMPO  ( 1963) 7 

100 

9.4UN 

MOD  (19  56-66)5 

0 

9.. ION 

Vestine  (1915) 17 

0 

9.40N 

NOO-OGS  1700  (19t,5)6 

0 

9.30N 

The  coincidence  of  the  geomagnetic  dip  equator  at  ionospheric  height 
and  earth’s  surface  lor  southeast  Asia  in  the  vicinity  of  longitude  100° 
east  reveals  that  the  earth's  magnetic  field  is  uniform  in  this  region. 
Thus,  the  dipole  field  model  is  a  good  approximation  of  the  earth’s  mag¬ 
netic  field  for  this  region. 

D.  Accuracy  of  Results 

A  constant  ionospheric  centroid  height  was  assumed  for  the  position 
of  the  transverse  propagation  in  the  ionosphere,  since  at  this  point  the 
magnetic  dip  angle  was  to  be  obtained.  In  reality,  the  ionospheric  cen¬ 
troid  height  is  not  a  constant  value.10  For  a  variation  of  centroid 
height  of  50  km,  errors  in  the  values  of  magnetic  dip  angle  and  the  cor¬ 
responding  suui onosphei ic  point  of  less  than  3  percent  will  be  introduced. 

In  this  experiment,  the  time  of  the  satellite  at  the  transverse  posi¬ 
tion  was  observed  on  analog  chart  records.'  The  accuracy  of  this  observa¬ 
tion  was  limited  by  the  width  of  the  recorder  marking  lines.  The  smallest 
division  of  time  on  the  marking  scale  was  5  seconds.  The  events  time 
record  can  be  resolved  to  within  1  second.  The  error  due  to  the  time 
reading  is  unlikely  to  exceed  2  percent. 

Since  the  magnet  declination  in  Thailand  is  small,5  in  the  complete 
Eq.  (A- 16 ) ,  this  term  is  neglected  for  Eq.  (11)  to  determine  the  magnetic 
dip  angle.  The  error  due  to  this  approximation  is  less  than  2  percent. 

However,  the  total  errors  in  the  values  of  magnetic  dip  angle  at  the 
altitude  of  ionospheric  centroid  height  obtained  using  the  beacon  satel¬ 
lite  technique  are  usually  less  than  10  percent. 
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VI  SUNMA  V  AND  CONCLUSION 


The  use  of  simple  instrumentation  for  Faraday  rotation  observations 
has  provided  a  technique  for  calculating  the  magnetic  dip  angle  o  er 
Thailand.  The  magnetic  dip  angle  at  ionospheric  heights  obtained  from 
the  Faraday  rotation  technique  is  very  close  to  the  surface  value  and 
compares  well  with  values  obtained  Sy  other  workers.  The  results  also 
have  been  compared  with  the  spherical  harmonic  function  with  good  agreement. 

This  technique  proved  useful  for  estimating  the  location  of  the 
geomagnetic  dip  equator  which  is  located  to  the  south  of  Bangkok.  It  is 
found  that  the  location  of  this  equator  at  ionospheric  height  is  over 
the  geographic  latitude  of  9.30N— about  480  km  south  of  Bangkok.  This 
value  is  very  close  to  the  surface  value  as  obtained  from  several  sources. 

Knowledge  of  the  coincidence  of  the  magnetic  dip  angle  at  ionospheric 
height  and  the  earth’s  surf-vc  i  '’cates  that  the  earth’s  magnetic  field 
in  Southeast  Asia  is  uniform.  ..*S  uniformity  of  the  earth’s  magnetic 
field  somewhat  resembles  that  of  a  uniformly  magnetized  sphere,  since 
this  would  conform  with  a  simple  dipole  field.  from  this  indication,  a 
dipole  field  model  is  a  good  first-order  approximation  of  the  earth’s 
magnetic  field  in  the  Southeast  Asia  region. 

The  results  of  the  work  presented  in  this  report  and  the  further 
analysis  of  these  satellite  data  for  their  ionospheric  data  will  give  a 
potential  of  improvement  in  understanding  of  the  radio  propagation  effects 
related  to  geomagnetic  control  of  the  ionosphere  in  Southeast  Asia.  Thus 
the  simple  instrumentation  used  for  Faraday  rotation  observations  can 
provide  part  of  the  ionospheric  data  base  required  for  the  study  of 
equatorial  radio  communication. 
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APPENDIX  A 

DERIVATION  OF  MAGNETIC  DIP  ANGLE  USING 
FARADAY  ROTATION  TECHNIQUE 
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APPENDIX  A 


DERIVATION  OF  MAGNETIC  DIP  ANGLE  USING 
FARADAY  ROTATION  TECHNIQUE 


Assume  no  retraction  in  the  ionosphere  or  lower  atmosphere.  Con¬ 
sider  the  geometry  shown  in  figs.  A- 1  and  A-2  where  the  coordinates  of 
points  A  (sub-satellite  point),  E  (sub-satellite  point),  and  ft  (ground 
station)  are  6 As;  6 f ,  A(1  and  t^r,  A r;  geographic  latitude  and  longitude, 
respectively.  The  earth’s  ladius  is  denoted  by  ft,  and  the  ionospheric 
centroid  height  is  h f .  Satellite  height  is  h s . 

By  using  spherical  trigonometry, 

cos  C  =  sin  ds  sin  6>r  +  cos  &s  cos  8r  cos  (As  -  A  r)  .  (A-l) 

Let  the  slant  distance  from  satellite  to  ground  station  be  equal 

to  l 


l  =  [ft,2  ♦  (ft,  ♦  ft,)2  -  2ft((ftt  +  K t)  cos  C]K 


( A- 2 ) 


The  inclination  oi  signal  ray,  at  the  ionospheric  centroid  height 

can  be  determined  as 


cot  I. 


fl,(ft,  ♦  hs)  sin  C 


{(Re  *  hc)2[R 2  +  (ft,  +  hs)2  -  2 flf(fte  ♦  hs)  cos  C ]  -  [Rf(Rt  +  ht)  sin  C]2}% 

(A-  3 ) 

From  the  ray- path  geometry  of  Fig.  A- 1 , 

I” _ flf(flf  *  ht)  sin  C _ "I 

[(fi,  ♦  h  c)  [ft2  ♦  (flf  ♦  hs)2  -  2  fte(ft,  ♦  A,)  cos  C]5*  J 


COS 


j  (ft,  *  li,!  sin  C  "j 

Jft2  +  (fle  +  A,)2  -  2ft,(fte  ♦  A,)  cos  C] '•  j 


.  (A- 4) 
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FIG.  A-2  SURFACE  ARC  DISTANCE 


Figure  A-2  shows  the  arcs  of  the  earth’s  surface  from  the  geograph¬ 
ical  north  pole  to  the  sub-satellite  point,  sub- i onospher i c  point,  and 
ground  station.  (If  more  convenient,  we  may  join  them  to  the  geographical 
south  pole.)  Then  the  sub- i  onospheric  point  i*(<?  ,A.  )  is 

b  (  =  sin”1  [sin  8  (cos  b  -  sin  b  cot  C)  +  sin  6 s  sir  6  cosec  Cl 

(A-  5) 

K t  =  Ar  +  N  (If  \f  >  Ar)  ( A - 6 ) 

»  \r  -  N  (If  Kt  K , )  ( A- 7 ) 

where 


/ 

>i"| 

■  -  1  d 

1  si  n  b 

1  I  — 

b s  -  sin  b r 

COS 

-V 

sin  ' 

— 1 

f 

COS  cf  t 

\ 

cos  b r  sin 

c 

/  _ 

w 

J 

G1 


The  declination  of  the  signal  ray,  Dg,  at  the  ionospheric  centroid 
height  can  be  determined  by 


D  K  *  cos'  1 


cos  6  sin  (C  -  6) 


-  tan  fcL  cot  (C  -  6) 


( A- 9 ) 


The  angle  between  signal  ray  path  and  earth’s  magnetic  field,  tf,  can 
be  determined  by  the  geometry  as  shown  in  Fig.  A-3.  Let  OX,  OF,  and  OZ 
be  three  perpendicular  lines  through  the  ionospheric  centroid  height  at 
point  0,  where  Jg  and  Dg  are  magnetic  dip  angle  and  declination,  respec¬ 
tively.  Using  the  direction  cosines  theorem,  18  the  arth’s  magnetic  field 
Ob  is  in  the  direction  with  the  direction  cosines  (cos  ag,  cos  / ig ,  cos  yg), 
and  the  signal  ray  OS  is  in  the  directions  (cos  ag,  cos  /ig,  cos  yg).  By 
the  rule  lor  the  addition  of  direction  cosines,  the  angle,  6,  between  OB 
and  OS  can  be  determined  by 


6  =  cos  ol  cos  ag  +  cos  cos  /3g  +  cos  y„  cos  v 


(A- 10) 


By  the  spherical  trigonometry,  we  can  solve  for  te  cosine  value  as 
follows: 
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cos 
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cos 

I  R 
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Dr 
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sin 

*r 

( A- 11) 

cos 

aB 
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cos 

*B 

cos 

Dr 

cos 

Pb 

= 

cos 

lR 

sin 

Dr 

cos 

yB 

= 

sin 

*B 

. 

(A- 12) 

Substituting  Eqs.  ( A- 1 1 )  and  (A- 12)  into  (A- 10) 
cos  6  =  cos  7g  cos  Jg  cos  (Dg  -  Dg)  +  sin  J 


R  sin  Jl 


(A- 13) 
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X  (  NORTHWARD  ) 


FIG.  A-3  RECTANGULAR  COORDINATE  OF  GEOMAGNETIC  FIELD  -- 
RAY  PATH 
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Note  that  I R,  lR  is  positive  when  the  direction  is  below  the  hori¬ 
zontal.  Since  Fig.  A-l  shows  that  I (direction  of  signal  ray  from 
satellite)  is  always  negative,  then  Eq.  (A- 14)  becomes 

cos  8  =  cos  I B  cos  /fi  cos  (Dfl  -  Dg)  -  sin  /fi  sin  / B  .  (A- 14) 

At  the  position  where  the  signal  ray  from  satellite  is  perpendicular 
to  the  earth’s  magnetic  field,  8  -  90°,  we  have 

cos  / cos  I g  cos  (Dfl  -  D fi)  -  sin  / sin  IB  =  0  ( A- 15) 

or 

I g  =  tui"1  [cot  Jg  cos  (Dg  -  Dg)]  .  ( A- 16) 

Since  the  magnetic  declination,  Dfit  is  small  in  Thailand  (always  less 
than  1  degree),5  the  approximation  can  be  written  as 

cos  (Dg  -  Dg)  =  cos  Dg  .  ( A- 1 7 ) 

Then  Eq.  ( A- 16)  can  be  rewritten  with  less  than  2  percent  error  as 

I fl  =  tan"1  (cot  / g  cos  D R)  .  ( A- 18) 

Equation  (A- 18)  is  used  to  calculate  the  magnetic  dip  angle  above 
the  earth’s  surface  from  the  selected  satellite  fading  records. 


APPENDIX  B 

FIRST  SIX  SPHERICAL  HARMONICS 
FOR  THE  GEOMAGNETIC  FIELD 


35 


HANK  PAGE 
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FIRST  SIX  SPHERICAL  HARMONICS 
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ABSTRACT 
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magnetic  field  vector  is  90°  (transverse  position). 

The  magnetic  dip  angle  calculation  was  based  on  the  assumption  of  a  constant  height  of  350  km 
for  the  centroid  of  the  ionospheric  electron-density  profile.  The  spherical  harmonic  analysis 
is  used  to  verify  the  experimental  results.  The  results  compare  well  with  the  surface  magnetic 
dip  angle,  measured  in  1956-1960  by  the  Geodesy  Department  of  the  Ministry  of  Defense,  Thailand, 
and  indicate  that  the  magnetic  dip  angle  at  ionospheric  height  i3  very  close  to  the  surface  value. 

These  data  permit  estimation  of  the  position  of  the  geomagnetic  dip  equator  at  the  ionospheric 
height.  The  estimated  position  of  this  equator  is  directly  above  the  geographic  latitude  of 
9.3fl°N.  This  latitude  intersects  the  southern  peninsula  of  Thailand  about  480  km  south  of 
Bangkok  (near  Surat  Thani).  This  estimated  value  compares  well  with  the  surface  geomagnetic 
dip  equator  and  indicates  tjjiat  the  dipole  field  is  a  good  model  for  Southeast  Asia  in  the 
vicinity  of  longitude  100°E. 
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